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Abstract. We present the results of phase-referenced VLBA+EfFelsberg observations at five frequencies of the double-image 
gravitational lens JVAS B02 18+357, made to establish the precise registration of the A and B lensed image positions. The 
motivation behind these observations is to investigate the anomalous variation of the image flux-density ratio (A/B) with 
frequency - this ratio changes by almost a factor of two over a frequency range from 1.65 GHz to 15.35 GHz. We investigate 
whether frequency dependent image positions, combined with a magnification gradient across the image field, could give rise 
to the anomaly. Our observations confirm the variation of image fl ux-density ratio with fr equency. The results from our phase- 
reference astrometry, taken together with the lens mass model of IWucknitz et alj j2004 . show that shifts of the image peaks 
and centroids are too small to account for the observed frequency-dependent ratio. 



1. Introduction 

The radio source JVAS B02 18+357 was identifi ed as a grav- 
itationally lensed system by IPatnaik et al.l ( Il993l) . The source 
(Figure consists of two compact images, A and B, sepa- 
rated by ~330 mas, the smallest angular separation amongst 
the known galaxy lenses, and a faint Einstein ring of a sim- 
ilar diameter. Optical observations of the lens (Brow ne et al.l 
Il993[ LStickel & Kuh3ll993[ lO'Deaet al.lll992l) have revealed 
molecular absorption lines associated with the lens galaxy 
at the redshift of 0.684. The morphology of the lens galaxy 
is spiral as suggested by the large differential Faraday rota- 
tion (~ 9 00 rad m"^) which is measured between the images 
JPatnaik e t al. 199 3 [) and detections of radio molecular ab sorp- 
tion lines dWiklind & Combeslll995t ICarilli et al.lll993h indi- 
cating that the lens-galaxy is rich in ionized gas, common to 
late-type spirals. The small image separation, which is indica- 
tive of a low-mass galaxy acting as a lens, lends support to 
this morphological categorization. The background source has 
prope rties of a powerful rad io source and a red optical contin- 
uum dStickel & Kuhi<ll993h . based on which it is conjectured 
to be a blazar, and its redshif 't has been confirrned u sing opti- 
cal sp ectroscopy (z - 0.944. Icohen et al1l2003t iBro wne et al. 
I1993I) . Further, it is variable in its radio emission and an accu- 
rate value for the time delay be tween varia tions in the images 
has bee n measured by Bi ggs et al.l ( Il999l) . (10.5 + 0.4) days 
(see also lCohen et al.l2000l) . 

Previous VLB I observations at high fr equencies 

llBi ggs et al.l l2003[ iKemball et all l200l[ IPatnaik et al...l995.) 
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Fig. 1. A combined EVN/MERLIN 1.7 GHz map of B0218+357, ob- 
served on 29'*^ May 1997 (Patnaik et al., unpublished). The reso- 
lution of the map is 50 mas and it shows an Einstein ring and the 
two compact images, A (right) and B (left). The contour levels are 
6.5 mjy beam"' x (-0.5, 0.5, 12, 4, 8, 16, 32, 64) and the peak flux- 
density is 650 mJy beam"'. 



show that both images consist of a compact core component 
and a secondary compact component, offset by ca. 1 .5 mas 
in the direction of a lower brightness jet of some 15 mas 
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Table 1. The image flux-density ratios as a function of frequency. Sp^ and Sg denote flux-densities of image A and B, respectively. Entries 
marked with an * indicate MERLIN and VLA measurements free from contamination by the Einstein ring. The error bars are included where 
available. 
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variability and the time delay; see text) 



extent. At lower frequencies (e.g. lPorcas & Patnaildll993) the 
resolution is insufficient to separately resolve the two compact 
components but the total extent of both images increases 
considerably with increasing wavelength. In particular, image 
A, the stronger and larger, is extended in PA -40°, the expected 
direction of "tangential stretching". 

The wealth of data coming from numerous radio and op- 
tical observations of this system at various frequencies and 
epochs provides constraints for lens models and, together 
with the time delay, give the possibility to determine the 
Hubble cons tant. Ho. This has been accomplished by vari- 
ous groups. iBiggs et all (11999') obtained a value of //o = 
69^[g km s"' Mpc"' u sing a singular isotherm al elliptical (SIE) 
mass distribution and IWucknitz et alJ ( l2004l) obtained a value 
of Hq = 78+^ km s^' Mpc^' using a SIE potential (SIEP). 
The only problem in Hq determination is related to the fact 
that the lens galaxy has weak optical emission, and is almost 
completely masked by the emission from image B (which is 
brighter than its counterpart in the optical wavebands). Thus, 
the position of the lens galaxy is poorly known relative to the 
lensed images and has been taken as a free parameter by the 
earlier authors in their (respective) lens models. A direct at- 
tempt to extract the optical centre of the lens-galaxy has been 
made by lYork et all (l2005l) who after subtracting the emis- 
sion from the images, A and B, describe the lens centre as a 
point about which the residuals are most symmetric. From their 
observations alone and using a SIEP lens model, they derive 
Hq = 66+^ km s ' Mpc ' or Ho = 19^1^ km s"' Mpc ' depend- 
ing upon whether the spiral arms of the galaxy are masked or 
not. 

Even though the basic lensing characteristics of this sys- 
tem are well-understood and well-reproduced by the current 
lens models, there are a few complications. One of them is the 
steady and systematic dechne in the ratio of radio image flux- 



density of the A and B images with decreasing frequency. This 
is in direct contradiction to the expected achromatic behaviour 
of gravitational lensing for point sources. Table 1 lists the ob- 
served ratios from previous observations of B02 18 -1-357 at var- 
ious frequencies and times, and using different interferometric 
instruments. The ratio (A/B) varies from 3.7 to 2.6 over the fre- 
quency range 15.35 GHz to 1.65 GHz. At frequencies where 
the source is variable, the ratio measured at a single epoch will 
be in error if the timescale for significant variations is compa- 
rable with the time delay. The three-mo nth VLA monitorin g of 
B0218+357 at 15 GHz and 8.4 GHz bv lBiggs et af](ll999l) has 
revealed variations in the image flux-densities on a timescale 
of (70-100) days. This, along with the measured time-delay 
between the images, can induce variations in the ratio at the 
level of (10-1 1) %. This is in good agreement with the 15 GHz 
ratios of the image flux-densities measured at the same epoch 
which fall in the range from 3.45 to 3.96, with the mean value 
of 3.74 and an rms of 0.02. At 8.4 GHz, the level of the ob- 
served variation in the ratio reduces to (7-8) % and the rms 
reduces to 0.01 . It is indeed expected that the variability ampli- 
tude of the background radio source reduces with frequency as 
the contribution from non-variable radio components becomes 
larger. Hence, although the time-variability in the background 
source emission along with the time-delay between the lensed 
images results in a spread of the estimated image flux-density 
ratios at a given frequency, it is clearly not the cause of the 
observed chromaticity in the flux-density ratio. 

We note that interferometric measurements may underes- 
timate the individual image flux densities if the beam size is 
much smaller than the extent of the image. If this occurs pref- 
erentially in one of the images it will affect the flux-density 
ratio; since image A is larger than image B, this effect might de- 
crease the A/B flux-density ratio. Conversely, insufficient res- 
olution may lead to contamination of the B image by emission 
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Fig. 2. The cur ves indicate con stant relative magnifications (A/B) for the best fitting lens model with the lens position at = 260 mas and 
>'o = 1 17.5 mas i Wucknitzl2002h . The filled circles represent the sub-components, 1 (west) and 2 (east), and the arrows represent the expected 
direction of the shift in the centroid of the brightness distribution of the images at lower frequencies, i.e. towards sub-component 2 which forms 
the base of the radio-jet. 



from the Einstein ring. We have included entries in Table 1 
from short baseline instruments (VLA, MERLIN; marked with 
an '*') only at frequencies where the resolution is sufficient to 
separate the ring emission from image B, to minimize this ef- 
fect. 

Although gravitational lens imaging of a point object is in- 
trinsically achromatic, frequency dependent variations in the 
image flux-density ratio of an extended source are possible, 
provided its structure changes as a function of frequency and 
extends over regions of different relative magnifications. It 
seems quite possible that this occurs in B0218+357 since both 
conditions exist. Firstly, the small image separation leads to 
relative large changes in magnification across the extent of the 
images. For simple isothermal ellipsoid models of the lens po- 
tential, a shift of 15 mas in the position of a point-source image 
can produce a change in relative magnification from 4 to 2.5 
(see Figure|2ji. 

Secondly, the frequency-dependence of the B0218-)-357 
image structures is strong, with previous VLBI observations 
showing a marked increase of image sizes with decreasing fre- 
quency. Furthermore, it is common for the radio spectra of one- 
sided AGN jets such as that seen in B02 18-1-357 to steepen with 
distance from the nucleus, providing a natural mechanism for 
producing a frequency-dependent position. The position of the 
radio peak at the base of the jet is also expected to move away 
from the nucleus at lower frequencies - the "core-shift". 

However, instrumental limitations prevent these effects 
from being seen easily. One is the change with frequency of 
the resolution available from VLBI observations. Another is 
the loss of absolute position information when phase self- 
calibration is used to make VLBI maps, which prevents ro- 
bust registration of maps at different frequencies. Frequency- 
dependent position shifts should, in general, show up as a 



change w ith frequency of the sepa ration between different im- 
ages (e.g. lPorcas & Patnaiklll99^ . However, such differential 
shifts are hard to measure in some cases, especially when the 
shift is in the same direction in both images. Furthermore, any 
"feature" (e.g. the core) used to define the positions of the dif- 
ferent images may, in fact, include different fractions of the 
emission in the background source, since the images have dif- 
ferent magnifications 

In this paper we present the results of phase-referenced 
VLBI observations of the lensed source B02 18-1-357, made to 
establish an unambiguous registration of the structures of the 
radio images at different frequencies. Although the expected 
core shifts are relatively small (mas), the positions of the cen- 
troids of the image flux-density distributions may change with 
frequency by larger amounts, and these are more representa- 
tive of the position of the emission. It is therefore important to 
quantify their changes with frequency when investigating the 
origin of the frequency-dependent image flux-density ratios. 

2. Inverse phase-referencing 

VLBI observations suffer from the corrupting influence of the 
troposphere, ionosphere and instrumental uncertainties on the 
interferometer visibility phase. It is well-known that the use of 
phase self-calibration to eliminate these errors results in loss 
of the geometrical phases which contain information regarding 
the source position relative to the antennas, resulting in the loss 

of any information on the source position in the sky. 

A lternatively, the technique of phase-referencing jAIej 
can be used, wherein observations of the target and a 
reference source are alternated frequently, allowing mapping of 
the target and determination of its position with respect to the 
reference. If the reference source is sufficiently compact that its 
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position is achromatic, it can be considered as an astrometric 
calibrator. Then, the phase-reference observations can be used 
to make a correct registration of maps of the target source at 
different frequencies, to study any frequency-dependent struc- 
ture. 

Although this technique is most often used for mapping 
faint radio sources, the target in our study, B0218+357, is 
sufficiently strong (~ 1 Jy) that we can use "inverse phase- 
referencing". Here, the target is used as the phase-reference for 
determining the corrupting phases, and phase-reference maps 
of a point-like, achromatic astrometric reference source are 
made. One advantage is that relatively faint sources can be 
used as astrometric references, permitting the choice of sources 
closer to the target, which minimizes telescope drive times 
and reduces any difference between the tropospheric and iono- 
spheric phase corruptions of the target and reference. Another 
advantage is that multiple astrometric reference sources can be 
used to guard against the possibility that any single one may 
have chromatic structure. 



3. Observations 

VLBI observations were made at five frequencies, 15.35 GHz, 
8.40 GHz, 4.96 GHz, 2.25 GHz and 1.65 GHz, using the 
NRAO Very Long Baseline Array (VLBA) and the Eff'elsberg 
radio telescope (Eb). Two 14 h source tracks were used (13''^ 
to 14* and 14* to 15* January, 2002) and for each track 
(19 h to 03 h UT for Eb, 19 h to 09 h UT for the VLBA) 
we switched between two receivers every 22 min. For the 
first track, observations were alternated between 4.96 GHz and 
2.25 GHz/8.40 GHz (S and X dual band) using 1 bit sampUng, 
and for the second track between 1.65 GHz and 15.35 GHz 
(without Eb) using 2 bit sampling. The data were recorded on 
four 8 MHz basebands at all the frequencies with the exception 
of eight basebands at 4.96 GHz. All observations were made in 
a single (LHC) circular polarization mode, except the dual S/X 
observations which were made in RHC. 

Apart from observing B0218+357, three compact reference 
sources were observed (see Table 2) along with a fringe finder, 
B0234-H285. Scans on B0218+357 of 1 min 20 s duration were 
interleaved by scans of 2 min duration on each of the calibra- 
tors in succession. The position reference sources were selected 
from the NRAO VLA Sky Survey (NVSS) catalogue. From an 
initial sample of candidate sources, these three were found to 
exhibit flat spectra and the most point-like struc ture on the ba - 
sis of EVN and MERLIN 5 GHz observations frorcas"200l'). 
They all lie within two degrees of B02 18-1-357 and are stronger 
than 25 mJy at both 1 .4 GHz and 5 GHz, with spectral indices 
flatter than -0.47 (where the flux density, S oc v"; v being the 
frequency and a the spectral index) . 

The data were correlated at the VLBA correlator, with an 
output averaging time of 2.1 s and a frequency resolution of 
0.5 MHz, and further processed in AIPS, the Astronomical 
Image Processing Software package provided by NRAO. 
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Fig. 4. Maps of 0215+364 at 2.3 GHz. a) Inverse phase-referenced 
map with ionospheric phase caUbrations applied to the data, b) 
Without ionospheric phase cahbrations. Both the maps are plotted 
with a 5 mas restoring beam and with identical contour levels at mul- 
tiples of 3 mJy. c) Hybrid map obtained by means of self-calibration. 
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Fig. 3. The top and bottom panels show 1.65 GHz, 2.25 GHz, 4.96 GHz, 8.40 GHz and 15.35 GHz maps of images A and B respectively, with 
frequencies increasing from left to right, top to bottom. The origin (0,0) is defined by the position of the peak intensity found in the initial run 
of fringe-fitting. This is the centre of A at all frequencies except for at 1.65 GHz at which the peak is in image B, entailing a reversal in the 
labeling of the axes. The residual noise in the maps is < 500 /jJy beam"'. The restoring beams used are (7, 5, 2, 1, 0.5) mas in the order of 
increasing frequency. The images clearly manifest all the earlier observed lensing-characteristics. Image A is tangentially stretched at a PA ~ 
-40° and image B seems to be tangentially compressed. At 8.4 GHz and 15.35 GHz the images reveal further substructure, sub-components 1 
and 2 with ~ 1.4 mas separation, representing the core-jet morphology of the background source. 
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Table 3. The image flux-density ratios from this data set. h denotes the resolution and cr is the root-mean-square noise in the maps. 
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Fig. 5. a) Image flux-densities obtained from these observations, b) Image flux-density ratio from these observations (with cross symbols and 
dashed error bars) and previous observations (with plus symbols and solid error bars, see Table 1). 



Table 2. The position references, their separations from B0218+357, 
the correlated flux density on the short EVN baselines (S-EVN) and 
on the long EVN baselines (L-EVN) at 5 GHz. The ratio L/S is a 
measure of source-compactness. 



Source name 


Separation 


S-EVN 


L-EVN 


L/S 




(degrees) 


(mJy) 


(mJy) 




NVSS B0210+366 


2.03 


127 


85 


0.66 


NVSS B0215+364 


1.03 


140 


86 


0.61 


NVSS B0222+369 


1.30 


151 


91 


0.61 



4. Data analysis and maps 

An essential step in the reduction of phase-reference data con- 
cerns the assumption of similar tropospheric and ionospheric 
phase errors for the reference and the target source. For VLBI 
observations the extra paths through the troposphere and iono- 
sphere are very different for different antennas, especially be- 
cause the source is seen at different elevations. The mean phase 
errors due to the troposphere, which scale as /t ', can be esti- 
mated approximately for each antenna, and are taken into ac- 
count in the model used for correlation. The mean phase er- 
rors due to the ionosphere, which scale as A, are highly unpre- 
dictable, however, and thereby prohibit the use of any model 
at the time of coiTelation. Since the errors become pronounced 
at long wavelengths, it is necessary to apply phase corrections 
from an ionospheric model after correlation. In applying the 



ionospheric and tropospheric models, it is only the mean phase 
offset that is removed for every antenna. The terms correspond- 
ing to small-scale temporal and spatial fluctuations still remain 
and it is hoped they are the same for the target and the reference 
source. For these observations, we used the AIPS task TECOR 
to apply an ionospheric model, produced by the Jet Propulsion 
Laboratory (JPL), and which gives an estimate of the total elec- 
tron content as a function of longitude, lat itude and time (for 
more details see IWalker & Chatterj eel 19991) . 

In the further analysis standard reduction procedures 
were used, applying amplitude calibration derived from tele- 
scope radiometry measurements and using fringe-fitting, self- 
calibration, CLEAN deconvolution procedures and mapping. 
As mentioned earlier, B0218+357 is sufficiently strong that 
self-calibration procedures could be applied at all frequencies. 
Hybrid maps were made by cleaning two sub-fields simultane- 
ously, one for image A and the other for B. After initial maps 
had been obtained the fringe-fitting was repeated using these 
as an input model, and the mapping was repeated. Full resolu- 
tion hybrid maps of both A and B images at all five frequencies 
are shown in Figure |3l 

All three reference sources also proved to be strong enough 
to use phase self-calibration, and hybrid maps were there- 
fore also made for these. On inspection, two of these sources 
were deemed to be unsuitable as position references. For 
B0210+366 the flux density drops to about 30 mJy at 1.6 GHz 
and consequently the position of the emission peak is not well 
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defined. The maps of B0222+369 reveal a jet like feature at 
2.3 GHz and 1.6 GHz, casting doubt as to whether the posi- 
tion of the peak is achromatic. Therefore for further analysis 
only B0215+364 was used as a position reference; it appeared 
point-like in the maps at all frequencies. An amplitude self- 
calibration procedure was used on this source to determine cor- 
rections to the a priori amplitude calibration, and these were 
then applied to the B02 18+357 data to improve the maps. 

Phase-reference maps of B0215+364 were made at all fre- 
quencies by applying the self-calibration solutions obtained 
from B0218+357 to the B0215+364 data. The dynamic range 
of these maps is poor and deteriorates with increasing fre- 
quency from 30:1 to 2:1 with rms noise ~ (1-1.5) mJy. 
To investigate the efficacy of the ionospheric model we com- 
pared phase-referenced images made with and without iono- 
spheric phase corrections. One such comparison can be seen in 
Figure El which shows maps of B0215+364 at 2.25 GHz. The 
peak flux per unit beam is almost 1 .5 times higher if an iono- 
spheric model is used [Figure |3 a)]. Moreover the peak posi- 
tions differ by almost half a milliarcsecond. We conclude that 
even though the rms deviation in the total electron content mea- 
sured by different groups is on the order of 25 %, the effect of 
applying any of the models is in the direction that improves the 
results. 

5. Results 

The image flux-densities and flux-density ratios from these ob- 
servations are presented in Table 3. To guard against loss of 
flux density due to over-resolution, we made low-resolution 
maps by discarding data from the longest baselines (using u,v 
resolution cut-offs). The flux densities from these maps were 
estimated by putting a box around the images and integrating 
the flux density within the box. The errors on the integrated 
flux-densities given in the table are derived from examining 
the spread in estimates obtained by varying a number of pa- 
rameters in the imaging process, and are larger than the formal 
errors obtainable from a Gaussian fit. Shown in Figure 13 a) 
are the image spectra and it can be seen that image A flux- 
density remains almost constant (to within + 25 mJy) at the 
upper four frequencies but drops suddenly by about 130 mJy 
(20 percent) at 1.65 GHz. In contrast, image B shows a grad- 
ual, monotonic increase with wavelength with no sharp drop. 
We note that the variation of the image flux-density ratio, as 
shown in Figure (Sjb), is similar in trend to that found for pre- 
vious observations, although the range is even larger (~ 4 to 
~ 2), with the value at 1 .65 GHz differing most from previous 
estimates. 

The result of the phase-referenced observations is shown 
in Figure |6| which shows the variation with frequency in the 
centroid position of the radio emission in the images with re- 
spect to B0215+364. To determine an image centroid position, 
the CLEAN components were convolved with a low-resolution 
restoring beam, and the peak in the resulting map was located 
using AIPS task MAXFIT, which fits a quadratic function to 
the brightest area of the map to determine the position ofi"- 
set with respect to the map centre. The same procedure was 
used to locate the offset from the map centre of the peak po- 



Fig. 6. The change in the centroid position relative to the position- 
reference B0215+364, at five frequencies in image A (top) and image 
B (bottom). The error bars at 1.65 GHz, 2.25 GHz and 4.96 GHz are 
dominated by the position uncertainty in the B0218+357 maps due to 
the use of relatively large beams, and at 8.4 GHz and 15.35 GHz due 
to the position uncertainty of B0215+364 caused by the poor dynamic 
range of the phase-referenced maps at these frequencies. 



sition of the calibrator, 0215+364, but using a full-resolution, 
phase-referenced map. The position of the B0218+357 image 
centroid with respect to 0215+364 (ignoring the constant po- 
sition difference used for correlation) is then given by the dif- 
ference of these off'sets (B0218+357 - 0215+364). The same 
procedure was used for image B. The upper panel of the fig- 
ure indicates a relative shift of ~ 5 mas in the centroid of 
image A between 15.35 GHz at 1.65 GHz, although the shift 
occurs only at 1.65 GHz. For image B, on the contrary, there 
is no shift detected in the centroid of the brightness distribu- 
tion. From the Singular Isothermal Elliptical Potential (SIEP) 
model (Figure |2ji, we find that a shift of 5 mas may cause the 
flux-density ratio to change as much as ~ 0.5, and looking at 
Table 3, this coiTesponds to ~ 12.5 % change in the observed 
flux-density ratio. However, this depends upon the direction 
in which the shift has occurred. From the phase-referenced re- 
sults for image A, we infer that this direction roughly coincides 
with the constant relative magnification contours (more or less 
tangential with respect to the lens centre) and this shift corre- 
sponds to < 6 % change in the observed ratio. Therefore we 
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Fig. 7. The A and B peak-to-peak separation (plus symbols) as a 
function of frequency. At 15.35 GHz, separations for both the com- 
ponents Al-Bl and A2-B2 are marked. The A and B centroid-to- 
centroid separation is also shown (in circles) but the error bars are 
omitted for the sake of clarity. 



draw the conclusion that the measured shift with frequency of 
either image centroid positions is not sufficient to account for 
the flux-density ratio anomaly for B0218-f357. 

An advantage of imaging gravitational lenses is that any 
frequency-dependent position in the source can be seen as a 
frequency-dependence of the separation vector of the lensed 
images, without requiring the use of an external reference 
source. As mentioned in Section 1, this method is insensi- 
tive in cases where the position shift is along the same direc- 
tion in both the images, but the accuracy of the measurement 
of the A-B separation is much higher than for that between 
B0218+357 and B0215+364. The peak-to-peak image separa- 
tions at different frequencies, measured with MAXFIT in full 
resolution maps, are shown in Figure Since the images at 
15.35 GHz are resolved into sub-components 1 and 2, the sep- 
arations of both Al-Bl and A2-B2 are indicated in the plot. 
We see that the separation at 15.35 GHz (Al-Bl), 8.4 GHz and 
4.96 GHz remains roughly the same, implying that the "peak" 



is centred around component 1 in both the images. At 2.25 GHz 
and 1.65 GHz, the separation increases in declination and de- 
creases in right ascension, in a direction opposite to component 
2. This is surprising as one might expect the separation at lower 
frequencies to gradually coincide with the A2-B2 separation, 
reflecting relatively more prominent emission from component 
2 at the base of the jet. Shown in the same figure as circles are 
the centroid-to-centroid image separations which are inconsi- 
tent with the peak-to-peak separations, although the behaviour 
of the latter is in concert with the image-structures seen at low 
frequencies (see below). 

At 1 .65 GHz we identify another component in image A, 
separated by ~ 10 mas from component 1 and ~ -33 ° P.A., 
shown in Figure |8(a)| The origin of this newly-identified fea- 
ture (hereon component 3) is of great interest and we have in- 
vestigated how this can be part of the bacJtground source. On 
applying the lens model bv lWucknitz et al.l (l2004h . we deduce 
that in image B it should be ~ 3 mas from component 1 as in- 
dicated in Figure |8(b)| Unfortunately the resolution at this fre- 
quency is not enough (~ 7 mas) to resolve this separation. This 
can be seen by placing a component at the expected position 
of component 3 in image B, of strength equal to the peak in- 
tensity of component 3 in image A [Figure |8(c)) divided by the 
relative magnification ratio at this position. On convolving this 
with the restoring beam, its effect can be assessed by compar- 
ing a slice made through it and the components 1 & 2 with the 
slice made through the same points in the observed image. As 
shown in Figure |8(d)| including another component in image B 
analogous to component 3 in image A has no observable effect 
other than changing the centroid of the brightness distribution 
by distance less than the positional uncertainity. Therefore, we 
cannot distinguish whether component 3 is a distinct feature in 
the background source imaged in A by the smooth macro po- 
tential of the lens, or whether it is caused by some other mech- 
anism. The shift in the centroid position of image A at 1 .6 GHz 
can be attributed to the existence of component 3. 



6. Discussion 

We have successfully used the technique of inverse phase- 
referencing to investigate the frequency-dependence of the 
emission from the images of B02 18-1-357; we believe this is the 
first time in which a gravitational lens has been used as a phase- 
reference. We have established the change in the centroid posi- 
tion of the image brightness distributions over five frequencies, 
and also investigated the change in the positions of the emission 
peaks using inter-image astrometry. The shift in the centroid in 
image A, which is significant (~ 5 mas) only at 1.65 GHz, is in 
a direction along which the relative magnification is predicted 
to be constant. In image B no significant frequency dependent 
shift is detected in the position of the centroid. A reasonable 
assumption is that the relative image magnification at the cen- 
troid positions derived from the model gives a good measure of 
the expected image flux-density ratio. Thus we conclude that 
the changing magnification gradient across the images is not 
the main cause of the anomalous change of image flux-density 
ratio with frequency. 
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Fig. 8. a) Second maximum (component 3) detected at 1.65 GHz in image A at a separation of ~ 10 mas from component 1 & 2 (The 
components 1 and 2 are not resolvable at these frequencies, hence what is marked is their superposition), b) The predicted position of component 
3 in image B at 1.65 GHz, marked as a black dot about ~ 3 mas from component 1 & 2. The resolution at this frequency is about 7 mas. c) A 
slice through the components 1 & 2 and 3 in image A. d) The predicted and observed slices through the components 1 & 2 and the expected 
position of component 3 in image B. 



As can be clearly seen in Figure at 1.65 GHz there is 
a large amount of low brightness emission that extends out 
to ~ + 30 mas. In comparison, at 15.35 GHz the emission is 
dominated by the compact sub-components with a separation 
of ~ 1.4 mas. At lower frequencies the (larger) images extend 
over regions where lens models do, indeed, predict significant 
changes in the relative magnification, and it may be insuffi- 
cient to simply consider magnifications at the centroid posi- 
tions. This is because low frequency emission from different 
regions of the background source is magnified by very differ- 
ent amounts. The resultant average magnification is thus the 
integral of the (background source) intensity-weighted magni- 
fication over the image area. 



We plan to extend our study by applying the model from 
LensClean and evaluating the ratio of the (background source) 
intensity-weighted magnifications of image A to B, explicitly 
taking the observed extension of the background source into 
account. We also want to extend our choice of lens model to 
one with a non-isothermal mas s-radius profile. Th e motiva- 
tion for this is the observation bvlBi ffgs etal]ll2003l) related to 
the jet stretching factor, which seems to be larger in image B 
than A, as seen when projected back to the source plane. But 
the deviation from non-isothermality is expected to be small 
iWuckni tz et al, 2004 ) since low-resolution VLA observations 
of the Einstein ring are fitted quite well with isothermal models. 

The identification of a distinct secondary maximum in im- 
age A at low frequencies (component 3), the sharp downturn 



10 



Rupal Mittal et al.: VLBI phase-reference observations of the gravitational lens JVAS B02 18+357 



in the spectrum of image A at 1.6 GHz (not visible in im- 
age B) and the small shift in the 1.6 GHz peak position in 
A are further results from these observations with no obvi- 
ous explanation in terms of the expected magnification gradi- 
ents across the images. It may, therefore, be necessary to con- 
sider more elaborate mechanisms, such as mass sub-structure 
(milli-lensing), free-free absorption or refractive scattering in 
the lens galaxy, to account for all the observed features in 
the B0218+357 images. Milli-lensing can produce significant 
frequency-dependent changes in the flux-density of one (or 
both) of the images provided the source size is comparable to 
the Einstein radius of the perturber and changes appreciably 
with frequency relative to it. The source size changes by a fac- 
tor 30 over the observed frequency range and its interaction 
with the caustics of the perturber can bring about frequency- 
depen dent changes in the total magnification JPobler & KeetonI 
l2005h . Non-gravitational effects such as free-free absorption 
or scattering will cause flux-density perturbations that fol- 
low an inverse frequency power-law. The effects produced 
by the above mechanisms on the image flux-density ratio in 
B0218+357 are a subject in their own right and are beyond the 
scope of the present article (iMitta]»2006t in preparation). 
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